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Abstract. Solutions of the wave equation in a space-time containing a thin cosmic string are examined in 
the context of non-linear generalised functions. Existence and uniqueness of solutions to the wave equation in 
the Colombeau algebra Q is established for a conical space-time and this solution is shown to be associated to 
a distributional solution. A concept of generalised hyperbolicity, based on test fields, can be defined for such 
singular space-times and it is shown that a conical space-time is S-hyperbolic. 



1. Introduction 

Weak singularities have for some time been used to model physically plausible scenarios such as thin 
cosmic strings, impulsive gravitational waves and shell crossing singularities. Such singularities typically 
admit a locally bounded metric and well behaved curvature scalars but none the less are still classified 
as singularities due to a low differentiable metric resulting from topological defects. 

One such example is the conical singularity, used to model thin cosmic strings, resulting from writing 
down Minkowski space-time in cylindrical coordinates 

ds 2 = -dt 2 + dr 2 + r 2 dcf> 2 + dz 2 

and identifying = with <f> = 2ttA for A < 1. 

The singular behaviour of this space time becomes apparent when one parallelly transports a frame 
around a closed curve around the axis A and finds that it undergoes a holonomy of the same magnitude 
as the deficit angle 27r(l — A). For non-singular space-times with a C 2 ~ regular metric it may be shown 
that such holonomies are trivial (See Wilson and Clarke (1996) for a proof); it is therefore the case that 
there is a problem with differentiability. In suitable Cartesian coordinates 

x = t, x 1 — r cos(4> / A) , x 2 = r sm(4> / A) , x 3 — z 

the metric may be written as 



ds 2 = -dt 2 + dz 2 + i(l + A 2 )^ 1 ) 2 + (dx 2 ) 2 ) 



which, although it is locally bounded, admits directional dependent limits as one approaches the axis. 

The low differentiability of the metric at the axis does not prevent the curvature being calculated as 
a distribution. One can not directly use conventional distribution theory techniques because this would 
involve the evaluation of ill defined products of distributions which are not well defined (Schwartz, 1954). 
Geroch and Traschen (1987) introduced a class of regular metrics for which the distributional curvature is 
guaranteed to be well defined and showed that such a space-time must have distributional curvature with 
support on a submanifold of at most co-dimension one, which is clearly not the case for the cone. Instead 
regularisation techniques may be used such as those employed by Balasin and Nachbagauer (1993) and 
especially those using the framework of Colombeau's non-linear generalised function theory (Colombeau, 
1984) such as Clarke et al., (1996). Using such techniques it may be shown that the non-zero components 
of the energy-momentum tensor density are 

T\{-gY' 2 = TM-gY' 2 = -2tt(1 - A)8^{x\x 2 ) 
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This form of the energy momentum tensor is what justifies the use of this kind of space-time to model a 
thin cosmic string. This suggests that the interpretation of the curvature as a distribution could be one 
such criterion for regarding these quasi-regular singularities as physically plausible. 

Another question relating to the physical importance of a singularity arises out of the question of 
to what extent the singularity disrupts the evolution of Einstein's equations and therefore global hy- 
perbolicity. For space-times with a C 2 ~ metric (which guarantees the existence of unique geodesies) 
this question is answered by requiring the space-time to be globally hyperbolic (Penrose, 1979). There 
are however a number of space-times with a lower differentiability for which, although they may violate 
cosmic censorship, there still may be well-posed initial-value problems for test fields. This led to the 
proposal of the concept of generalised hyperbolicity (Clarke, 1998) in which one examined the extent 
to which singularities were obstructions to the evolution of the wave equation. In order to apply this 
concept to space-times with shell crossing singularities Clarke replaced the initial value problem for the 
wave equation 

Ou = f, u\ S = v, (n a d a u)\ s = w 

by a distributional version obtained by multiplying by a test field tp and integrating by parts once (rather 
than twice as is usual) to give 

/ d a ud b ipg ab fj,= ipffi- tpwfis 

JM+ Jm+ Js 

He then said that a space-time was □-globally hyperbolic if the above equation had a unique solution 
for all ip G T>(M). Moreover it was shown that this form of generalised hyperbolicity was satisfied for 
a class of curve-integrable space-times and in particular it was demonstrated for the shell-crossing dust 
space-times. 

In this paper we shall consider the question of generalised hyperbolicity for the conical space-time (1). 
Unfortunately this space-time does not admit a locally square integrable connection, so we are unable 
to apply the result of Clarke (1998) directly. Instead we shall follow a different approach: We shall use 
Colombeau's generalised functions to overcome the ambiguities which arise when attempting to multiply 
distributions which would arise when considering solutions to the wave equation in a space-time of such 
low differentiability. This involves writing down the Cauchy problem for wave equation in the Colombeau 
algebra Q (M) , and proving the existence of a unique generalised function solution. This will be done for 
a class of singular metrics which includes the 4-cone (1) (Section 4). Having obtained a unique solution 
we shall examine to what extent it is possible to interpret it as a distribution. We say that a space-time 
is (/-globally hyperbolic if there exists a unique solution to the wave equation in the Colombeau algebra 
Q and that this solution is associated to a distribution. We therefore demonstrate in this paper that a 
conical space-time is (/-hyperbolic. 

2. Colombeau's generalised functions 

We first briefly recall the essential details of Colombeau's theory (see e.g. Colombeau, 1984 or Ober- 
guggenberger, 1992 for further details). We denote the space of smooth functions with compact support 
2?(R n ), so that the space of distributions is T>'(M. n ) and we define the following spaces; 

Definition 1 . We define the space of smoothing kernels A q (MP ) as the space of all functions ipW G T>(MP) 
such that 

JR™ 

f ip (n) {0t k d n t = Vfc G N™ such that 1 sC |fc| q 

JR" 

where we are using a multi-index notation so that £ fc = (£ 1 ) fel . . . (£") n . 

Given a function / G L^^W 1 ), we may smooth it by constructing the following convolution; 

f(v>e,x)= I f(Q<pi n Hz-x)<rt 

where 

¥£°(0=e-V n) (£/e) 
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and for more general distributions T G V (R n ) we have 

f(<p e ,x) = (T,T x tpW) 

where 

We would like to construct a differential algebra of such smoothings; We define a base space £(R n ) as 
the space of all functions / : Aq x M™ — > C which are smooth as a function of x. Certainly smoothings 
of distributions are elements of this space; however one desirable property would be for the smoothing of 
/ G C°° to coincide with its natural embedding in that 

and 

f(ip e ,x) = f(x) 

are equivalent. This is achieved by constructing a space of null functions J\f(R n ) such that / — / G Af(R n ) 
and by working with a quotient space £(R n )/Af(R n ). 

Definition 2. The space of null functions Af(R n ) consists of the functions / G £ (R n ) such that given 
K C C R n and any collection of indices pi, ■ ■ ■ ,Pk with fc G N there is some N G N and an increasing 
unbounded sequence (7,) such that for each ip G A q for q ^ TV, 3c, 77 > such that 

sup |<9 P1 . . . d p J{<p e , x)\ ce 7 ^ (0 < e < rj) 

x£K 



The quotient space £(R n )/Af(R n ) is not well defined because 7V(K") is not an ideal of £(R n ). We may 
multiply a null function with a function of non-polynomial growth in to give a function which is not 
null. This problem is rectified by restricting £(R") to a space of moderate functions £m(R") which itself 
is an algebra and of which 7V(K") is an ideal. It may be verified that the smoothings of distributions are 
indeed moderate. 

Definition 3. The space of moderate functions £m(R") consist of the functions / G £(R n ) such that 
given K C C K™ and any collection of indices pi , . . . , pk with k G N there is some N G N such that for 
each Lp E An, 3c, 77 > such that 

sup \d Pl . . . d p J(tp E , x)\ ^ cs~ N (0 < e < i]) 



We finally define our space of generalised functions as the quotient space 



The space Q(W l ) is a differential algebra, of which C°°(R n ) is a subalgebra and X>'(R") is a linear subspace. 
In practice one usually works with the representative moderate functions to perform calculations in this 
space. Many elements of Q(R n ) do not result as the smoothing of distributions, but even so are equivalent 
to a distribution in the following sense; 

Definition 4. We say that [/] G Q(R n ) has an associated distribution T G V'(R n ) (written as [/] « T) 
if for any representative / and ip G T>(R n ) there exists some N G N such that 

lim [ f(ip e ,x)iP(x)d n x= (T,V) 



It may be shown, for example, that if H and 5 are the smoothed Hcaviside and delta distributions 
respectively then [H] 2 w H and [H][S] w \8. On the other hand, the generalised function [5] 2 has no 
associated distribution. 
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When one applies the theory to initial value problems one usually makes the additional requirement 
that the n-dimensional smoothing kernels are Cartesian products of one dimensional kernels (Oberguggen- 
berger 1992) so that 

V (n \x\... ,x n )=<pV>(x 1 )---<pV>(x n ). 

The reason for doing this is that it gives a natural way of defining Cartesian products and restrictions. 
However in order to simplify the calculations in Section 5 it is useful to work with rotationally symmetric 
smoothing kernels. Unfortunately it is not possible to write such a kernel as a finite sum of Cartesian 
products and for this reason we describe in Appendix A an alternative way of restricting 4-dimensional 
smoothing kernels to an initial surface S so that any spatial symmetries are preserved. Note however 
that the results of the next two sections are not sensitive to the precise restriction process adopted. 

3. The wave equation in the Colombeau algebra Q(M) 

We first introduce the convention that lower case Latin indices a, b, ... run over ... 3, Greek indices 
a, (3, ... run over 1 ... 3 and upper case indices A, B, ... run over 1 and 2. Suppose we have a space- 
time, equipped with a locally bounded singular metric, (M, g) and we wish to solve the Cauchy problem 
for the wave equation 

Du(t,x a ) = 

u(0,x a ) = v(x a ) 
d t u(0,x a ) = w(x a ) 

where initial data (v,w) lying in the Sobolev spaces 7i 1 (5) x Ti°(S) is prescribed on the initial surface 
5, t — 0. If a solution exists one would expect it to be defined as a distribution. This however will cause 
difficulty in interpreting 

Du=(-g)- 1 / 2 d a ((-g) 1 / 2 g ab d b u) 

as a distribution in the framework of Schwartzian distribution theory because it contains products of u 
with a singular metric and its weak derivatives which are not well defined in this theory The way to 
overcome these shortcomings is to employ the non- linear generalised function theory of Colombeau (1984), 
which does give a distributional interpretation to many distributional products that would otherwise be 
undefined. We first embed the metric g a b into the generalised function space Q{M) by constructing its 
representative [g e ab ) £ £m(M) as the convolution integral 

9l b {t, x a ) = J g ab (t + eC, x a + eC)^ 4) (C, O dC d 3 £ 

Since the initial data (v, w) does not have to be smooth, we must also embed it into the space Q{S) x Q{S) 
as (V, W) by defining the convolution integrals v £ and w £ in £m(S); 

Ve(x a ) 

w £ (x a ) 

The Generalised function wave operator □ : Q{M) — ► Q(M) then may be defined by denoting OU for 
U G G(M) as the class represented by (D £ u E ) £ £m{M) where (u s ) £ £m(M) is a representative for U 
and 

U'u e = (-g^dad-g^gfdbu,) 
We would like to then be able to solve the Cauchy problem in the space Q{M) 

DU(t,x a ) = 

U(0,x a ) = V(x a ) (2) 
d t U(0,x a ) = W{x a ) 

and obtain a solution U £ Q{M) which is associated to a distribution. In practice one works with the 
equivalent problem in £m(M); 

n^u e (t,x a ) = j £ {t,x a ) 

u E (0,x a )=v E (x a ) (3) 
d t u E (0,x a ) =w E (x a ) 



J v(x a +eC)^ 3 \C)d^ 
J w{x a +e£ a )<P {3) (C)d\ 
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where (f £ ) G Af{M) and obtain the associated distribution (if there is one) u by evaluating the limit for 
each ip e V(M) 

(u^) = \imj u e (<p,(t,x a ))i;(t,x a )» s (t,x a ) (4) 

where ip G A q for q G N large enough. 

The question of existence and uniqueness of solutions to (2) is not as straight-forward as for its 
classical function counterpart. Although one may establish the existence of unique solutions for (3) (for 
each ip G Aq an d e > 0) by using the using the Cauchy-Kowalewska theorem for the case of analytic data, 
this is not enough to establish the existence of unique solutions to (2). One must carefully examine the 
role that moderate and null functions have to play. (See Oberguggenberger (1989) for examples of similar 
Cauchy problems that do not admit unique solutions). Firstly, we want U to be an clement of Q(M). 
This means that we have to verify that (u £ ) G Em{M). Secondly for U to be a unique we only require 
that the solution (u £ ) of (3) arc unique up to an element of M{M). This means that for a particular 
choice of (f £ ), (v £ ) and (w e ) we are allowed to have more than one solution provided they are all well 
defined elements of £m(M) and differ by elements of J\f(M). We must also take into account that (f e ), 
(v £ ) and (w £ ) could have been chosen with freedom up to a null function, so we require that the solution 
(u £ ) is augmented by at most a null function if (/ e ), (v £ ) and (w £ ) are augmented by null functions. Since 
this Cauchy problem is linear it follows that the solution U G Q(M) is unique if for the problem (3), with 
(v £ ), (w £ ) G Af{S) implies that (u £ ) G M{M). 

Moreover even if there is a unique solution U G Q(M), it need not admit an associated distribution 
because the limit (4) may be dependent on the smoothing kernel if G A q or may not even be defined, 
however large we may choose q. 

4. Proof of the existence of unique solutions in Q(M) 

We start by observing that we may take a representation in £m{M) of the embedded metric which 
admits the form; 

ds e 2 = -dt 2 + dz 2 + g AB (x c ) dx A dx B (5) 
where g AB G £m(M) and locally g AB and its derivatives may be bounded as follows; 



\9ab\<M 

\g? B \^M 

|0Ci ■■■dc k g AB \ < 



with Mfe being positive constants independent of e. It is easily seen that on embedding the conical 
metric (1) into Q(M), the g AB and its derivatives admit this form; 

gl 1 (x A ) = \(l + A 2 ) + \(l-A 2 )P £ (x A ) 

gh(x A ) = l(l-A 2 )Q £ (x A ) 

gUx A ) = \(l + A 2 )-l(l-A 2 )P £ (x A ) 

where 

and if^ is a two dimensional smoothing kernel obtained by integrating out the t and z dependence. 

Our aim is to estimate solutions u £ of (3) and its derivatives in terms of powers of e, given the 
moderate and null bounds of f £ , v £ and w £ using a method of energy estimates following Hawking and 
Ellis (1973) and Clarke (1998). We shall assume that K is a compact region of M which intersects the 
initial hypersurface S = So described by t = and that and {S' t }o^t^ti is a family of smooth space-like 
hypersurfaces in K which intersect So on a common 2-surface. See Figure 1. 
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We shall assume that the surfaces S T may be expressed by the equation a(t, x a ) = r; we may therefore 
define a normal to S T , n a = d a T and a surface element fi s s which satisfies 

fi E = {-g E ) 1/2 d 4 x = h% t A<1t 
We shall write the region bounded by So and S T as 

C€[0,t] 

Since we will have to estimate u e and its derivatives of arbitrary order, we shall construct a hierarchy 
of energies. It will turn out that it is much more convenient to work with energy integrals based on 
covariant derivatives rather than classical Sobolev norms based on partial differentiation. We define our 
energy integrals by the formula 



.5=0 Jb ^ 



E £ fc (r) = > / T" Ej t b n a »% T (6) 



where 



T ab — (n ac n bd — Ln ab n cd \pP 1 1 1 pPk-lQk-l \7etr7s vy 6 11 X7 £ X7 e V7 e 11 

1 s,k — \9e 9s 29e 9e ) e e ■ ■ ■ e s V c v Pl ■ ■ ■ V p fc _i u £ V d V qi ■ ■ ■ V q k - t u e 



rpab _ _l n ab 2 
1 e,0 — 29s a e 



e,0 — 2»£ " £ 
*a = d a t 
ef = gf + 2t a t b 

The form of the metric (5) guarantees that V„t(, = This enables the tensor e" b to be constructed so 
that it is positive definite and is annihilated by covariant differentiation i.e. V^e bc = 0. 
We also define the following 3-dimensional Sobolev norms 

1/2 

• - Is. ~ "' "' ''' " 2 "" 



{ ^ / |a pl ...a w « e | 2 /iS T | 

^ Pl-.-Pj J S T ) 

Mt = { E / Iv^.-.v^lVsl 



1/2 



■ P1---PJ • 

These are equivalent to the energy integrals in the following sense 
Lemma 1. There exist positive constants A k , A' k , B k and B' k such that 

E k £ (r)^A k (\\u E \\ k s f (7) 
(\\us\\ k s f < A' k E k £ (r) (8) 

k 

(|| Ue |||J 2 ^B fe ^^i^(IMi) 2 (10) 

||it e ||s , ||u e |ls T an( i E*(t) are functions of t and We now define the concepts of having a moderate 
or null bound for such an object. 
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Definition 5. We say that p : Aq — > C has a moderate bound if there is some N G N such that for each 
<p 6 An, 3c, 7] > such that 

\p{ip e )\ ^ ce~ n (0<e<ry) 

Definition 6. We say that p : Aq — > C has a null bound if there is some JV € N and an increasing 
unbounded sequence (7 9 ) such that for each (/? £ for q ^ N, 3c, i] > such that 

|pfo> e )| < ce"-" (0<£<7?) 



An important consequence of Lemma 1 is that E^(t) has a moderate bound if and only if ||u e || <y has 
a moderate bound and that E*(r) has a null bound if and only if ||«e||| has a null bound. It is therefore 
sufficient to estimate these energy integrals. The idea is to form an energy inequality to give a bound on 
E^(r) in terms of the lower order energies El(r),... .E^-^t) and E*(0),. .. ,E*(0), which is determined 
by the initial data, and to show that the properties of the E^(0) having a moderate (null) bound may be 
carried through to E^(r). 

We first integrate around the region Q T and apply Stokes' theorem to obtain 

E e fc (T)-E^(0)=]T / t b V%T%p* 

3 =0 Jn r 

where 

VIT*\ = . . . e»-i«-i g^gf V^V^ . . . V^_ lUe V|V^ . . . V% k _^ e 

+ 2e™ . . . e^- 1 "- 1 5 f 5 f V|V^ . . . % k _u £ Vf a V £ d] . . . V^ £ 
V|T» 6 = -V* Ue « e (11) 

This expression involves derivatives of order k + 1, but we may eliminate these substituting the wave 
equation into it. For the base energy case k = 1; 

V £ a T« b = - We V^ £ 

Therefore 

E*(t)=E*(0) + / t a V> £ (/e-«e)M £ 

On making estimates; 

E^(r)^Ej(0) + iiKlJ, T (||/ e ^ T + rt T ) 

^E^w+iLidi/^jVi^dKi^) 2 

where 




1/2 



On applying (8) this inequality may be written as 

El{r)^^) + \Li{\\feil T ) 2 + -2 L ^ [KiQdC (12) 
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An energy inequality may be also obtained for higher order energies; we must first eliminate the 
derivatives of order k + 1 from the expression (11) by using the differentiated wave equation 



sfv* 1 ...v; fe _ 1 v*v^ = v* 1 ...v^_ 1 / £ 

In order to substitute this into (11) we must shuffle the covariant derivative indices of u by repeatedly 
applying the Ricci identities 

2^[a^fc] Xp 1 p 2 ...p k — Rabpi -^-cp2-..pk Rabp2 ^p\c...p k ~\~ ' ' ' ~\~ Rabp k ^p\p2...c 

In this way we find that that 

9 f vsvf . . . % k _u £ = v; x . . . % k _j e + E nf-}j}_ iUe (is) 

3=1 

2Vf a n, . . . V; k _ lUe = 11%^} u s (14) 

where TZ^'^Ug represents a linear combination of contractions of the (k — j)th covariant derivative of the 
Riemann tensor with the jth covariant derivative of u £ . This quantity may be bounded as 



£ 2(2+fe-j) 

q\...qj 



where Ck,j is a uniform constant. Now the expressions (13) and (14) may be bounded as follows; 

\ g f vmvi . . . v^i 2 < k\v pi . . . v^m 2 + kj2 Ht'dl* 

3 = 1 

^c k E |v^...v^ i/£ | 2 

qi---qk 



F 2(l+k-j) I 91 
91 — 9j 
l^j^fc-l 



|2Vf a V^V^...V^ iU£ | 2 <§ E |V^...V^ £ | 2 



<?l ■••<?fc-l 

where Ck is a uniform constant. Therefore 

1/2 



\^ a T%\^cd E l v pi--- v p^l 2 } 

Pl-Pfc 

r i i 1/2 

x E l^ 1 ---V^_ 1 / £ | 2 + E ^(TT^)I V pi---^^| 2 
l pi---Pfc-i Pi-Pi J 

e iv^...v^ £ i 2 + e iv^...v^ l/£ i 2 + e ^t^)I v pi---^^i 2 } 

'"Pl-Pfc Pi— Pic-l Pi— Pj ' 

l^j^fc-l 

On integrating this becomes 

E*(r) < E*(0) + cd(\\uMS + (ll/efc 1 ) 2 + E (IM^) 2 } 
which may be turned into an energy inequality by applying Lemma 1 

E*(r) < E*(0) + C'Shit 1 ) 2 + C 'k f E e(0 d( 

+ [mode (is) 

3=1 C "° 

Having obtained this energy inequality we shall prove the following 
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Lemma 2. Suppose that ||/ £ ||q , • ■ • , ||/e||q _1 have a null bounds. Then 

(1) if E°(0), . . . ,E^~ 1 (0) have moderate bounds then so does E*(r) 

(2) if E°(0), . . . , E£ _1 (0) have null bounds then so does E*(r) 

Proof. For part 1, we proceed by induction (Proof of part 2 is similar). For the case k = 1 we may apply 
Gronwall's Lemma to the energy inequality (12) 

El(r) ^ (El(0) + l Ll (\\f e f n y)e 3L ^ 

therefore if for ip G A q (q ^ N) it is the case that E*(0) = 0(e- N ) and ||/ e jj& T = 0(e 7 '- JV ) it will also 
be the case that E*(t) = 0(s~ N ) 

Now suppose we have for ip E A q (q ^ N); 

\\f s \\k=0(e^- N ) j = 0...k-l 

W £ (t) = 0(s- n ) 
E £ fe (0) = 0{e~ N ) 

on applying Gronwall's inequality to (15) gives 

I j =1 6 Jc=o ) 

therefore 

E £ fe (r) =0( e -( JV + 2fe )) 

We now want to apply Lemma 2 to show that if v e , w £ G £m{S) then u e G £m{M) and similarly if 
v £ , w £ E JV(S) then u e G Af(M). 

The initial data (u e ,ty e ) may be used to determine the initial energies E£(0) by differentiating the 
Cauchy problem (3) and using the wave equation as follows 

dpi--- d Pk u £ {0, x a ) =d pi ... d Pk v £ {x a ) 
d t d Pl . . . d Pk u £ (0, x a ) =d pi ... d Pk w £ (x a ) 

8ld pi . ..d Pk u £ (0,x a ) = 8 P1 . ..d Pk {gfd a d p diu e +g?TZ p d 7 diu e }(p,x a ) 

-di- 2 d Pl ...d Pk f £ (0,x a ) (16) 

and then using (6) to define the energy integrals 
Equation (16) may be expressed as 

k+2 

d{d Pl . . . d Pk u £ (0, x a ) = Yl G C1..2 (0. xa ) d °i ■ ■ ■ d°rt~ 2 Ue% x a ) 
i=i 

-d{- 2 d Pl ...d Pk f e (o,x a ) 

where Gf p \ is constructed by sums and products of the metric and its derivatives, hence is an element 
of £m(S). It therefore follows (by an inductive argument) that 

(1) v £ , w £ E £ M (S) and f £ E JV(M) implies that d{d pi . . . d Pk u £ (0, x a ) G £ M (S) 

(2) v £ , w £ E J\f{S) and f £ E N{M) implies that d{d pi . . . <9 Pfc u e (0, x a ) E Af(S) 

We therefore have proved the following 
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Lemma 3. Suppose that (v e ), (w £ ) G £m{S) and (/ e ) G M{S) then E^(0) /ms a moderate bound, and if 
in addition that (v £ ), (w £ ) G N(S) then E^(0) has a null bound. 

We also want to translate the bounds for E k (r), as given by Lemma 2 back to bounds for u e and its 
derivatives. This may be done by applying Lemma 1 in conjunction with the Sobolev embedding theorem 
(Hawking and Ellis, 1973). Suppose M is a manifold admitting an n dimensional embedded submanifold 
M' then according to the theorem given k > 0, 3M k > such that Vu G H k+m {M') with (2m > n) 

\d pi ...d Pk u\^M k \\4 k + m (17) 

where ^ 

lA = ( E / |5 P1 ...9 P3 «|V4 (18) 

< P!...pjJ M 1 ) 

and H k (M') is the Sobolev space of functions it existing almost everywhere and for which || 

It should be carefully observed that this theorem only gives bounds on the derivatives in the tangential 

directions to M' and likewise it shows that they may be bounded above by a Sobolev norm that only 

involves those derivatives. 

We first of all apply (17) to u £ on the submanifold S T . This gives us for u £ 

\d Pl ...d Pk u e \ ^ PoMK\\ k st m ^ P°M\us\\ k st m (2m > 3) 

Where the second inequality comes from the fact that we may replace the tangential Sobolev norm with 
the more crude version as its upper bound. However the above result does not put any bounds on 

any of the time derivatives. We do this by applying (17) to d(u £ ; 

\d Pl ...d Pk dlu £ \^P,. k \\dlu £ \\ k + m 

Therefore we have shown that 3M k > such that 

\d Pl ...d Pk u £ \^M k \\u £ \\ k + m , (2m>3) (19) 

On combining (19) with Lemma 1 we have proved the following 
Lemma 4. 

(1) IfE k (r) has a moderate bound then u £ G £m(M) 

(2) 7/E^(r) has a null bound then u £ G W(M) 

On combining Lemmas 2, 3 and 4 we have proved our main theorem. 
Theorem 5. 

(1) Ifv £ , w £ G £ M (S) and f £ G M{M) then u £ G £ M (M) 

(2) Ifv e , w £ G M(S) and f £ G N{M) then u £ G Af{M) 

We therefore conclude that for metrics of the form (5), a unique solution U G Q{M) exists to the 
Cauchy problem (2). 

5. The distributional interpretation of the solution 

It has been established that the solution u £ may be interpreted as a representative for an element of 
Q(M). We would like to now discuss whether or not it may be interpreted as a distribution. Thus we 
want u £ to be associated to a distribution u; that is for each ip G T>(M), 3q G N such that the limit 

(u,V) = lim / u £ (ip,x)ip(x)n £ (x) (20) 

is well defined and independent of ip. 
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To do this we consider a representation for u £ in terms of a Greens function. Assuming that u £ = 
for t < 0, we may write the Cauchy problem in an 'inhomogeneous' form 

□X + = It 

where 

u+(t,x a )=H(t)u £ (x a ) 

ft(t,x a ) = H{t)f £ (x a ) - 5'(t)v £ {x a ) - 8{t)w £ (x a ) 
In this way one may express the solution in terms of a Greens function (see Appendix B) 

<(x) = J GtM).ft(0v E 

It is therefore the case that if the limit (20) exists it will be given by 

(w,V>) = lim / u+(ip, x)il){x)if (21) 



M 



On expressing uf in terms of the Greens function and using the self adjoint property of □ together 
with the fact that G + is integrable on M x M, the order of integration may be interchanged giving 



M 



u+{ ! p,x)1>{x) ll = [ ( [ Gt. x (0ft(0m)^h £ (x) 
Jm \Jm J 



J MxM 



G-{x^)ft{x)^)fjf{x,0 

MxM 



I 

J A 



\ £ (x)f £ (x)n e (x) 

M 

- f (d t X £ (0,x a )v £ (x a ) + X £ (0,x a )w £ (x a ))^ s (x a ) (22) 



where 



k(x) = I G- x (mo^(o 

Jm 



Since f £ is null the first term on the right hand side of (22) will vanish as e — > for large enough q G N, 
provided that A e admits at most a moderate growth. For the second term to admit a well defined limit, 
independent of f e A qi as e — > 0, we require that the limiting function of X £ is also well defined and 
independent of ip. This is because v £ and w £ are the embeddings of locally square integrable functions 
on S. Thus to show that u £ is associated to a distribution we need to show that A e is sufficiently well 
behaved. To do this we follow the approach of Bruhat (1962) and use the fact that for e > 0, X £ (x) is a 
solution of a Volterra type integral equation involving the biscalar K (see Appendix B for details) 



K{x) 



±- ( UZK £ {^x)X £ {i)^(i) = ±- f KefoxMQMt) (23) 

27T J C -(x) ^ JC-{x) 



where /Xp is the volume element induced by \i e on C~(x). We next examine the limiting behaviour of 
the terms in (23). We first consider the integral on the left hand side. In Appendix B it is shown that in 
normal coordinates based at £ 

LTif £ (£, x) = -\K e % x)R e (0 + 0(e) (24) 
On the other hand using the results of Clarke et al (1996) 

iim| R £ (o*(om = mi - A ) j *(e, o,o, e)dfde (25) 
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So that provided Ao is well defined then 

lim / \yK e (£,x)\ e (t)i4(£) = -l(l-A)ir [ Aq(Oma-(0 (26) 

where A~(x) is the limit of the intersection of C~(x) with the axis A and fi^- is the volume form induced 
on A - (a;). It is shown in Appendix C that these quantities are well defined. 

The integral of the right hand side of (23) is more straightforward since for x, £ ^ A, K E (£,x) — > 1 
as e — + 0. Also we show in Appendix C that C~ (x) tends to a well defined limit Cq(x). However some 
care must be taken in interpreting Cq(x). One cannot simply take this to be the past null cone of x in 
the conical space-time with the axis removed. This is because deleting the axis results in a tear in the 
null cone and destroys the S 2 topology. However (in the case of a rotationally symmetric smoothing) 
it is shown in Appendix C that Cq (x) — lim e ^o C~ (x) is well defined and does have S 2 topology the 
missing piece being generated by geodesies which pass within 0(e) of the axis A, and which hit the axis 
in the limit. One can also show that the corresponding volume form p^(0 also has a well defined limit 
for £ ^ A given by the volume form induced by the conical metric on Cq (x) . Furthermore the value on 
the axis is bounded so that we need not include the contribution from the integral over the axis. Hence 
provided that Aq exists we have 



lim 



/ K s (t,x)m»uo= [ ^(0/4(0 (27) 

*°Jcr(x) Jc^(x) 

where Cq(x) = Cq(x) \ A. 

On the basis of (26) and (27) we now define X(x), to be the solution of 

A(s) - ^(1 - A) f Ao(0 MO = ^ / m /£(0 (28) 

Note that this gives A in terms of ip on Cq (x) and A on Cq (x) n A. For a point x that lies on the 
axis A~ (x) degenerates into a pair of null lines and /za vanishes so that the integral on the left hand side 
of (28) vanishes and we have 

\(x) = ^-[ V(0Mr(0, *eA (29) 
2tt Jc-( X ) 

where for a point x on the axis, the past null cone Cq (x) is given in quasi-Cartesian coordinates by the 
usual Minkowskian formula. Equation (29) gives an expression for A G A simply in terms of an integral 
involving i[>, so we may substitute back for A in the left hand integral of (28) to obtain an integral for A 
in terms of ip which is valid for all x. 

A(z) = ^-/ m4(0 + ^-(i-A) [ ([ m» m w)»A(t;) (so) 

27T J£ec-(x) D7T JzeA-(x)\J v ec-(i) J 

Thus a solution to (28) is given by the above integral (30). Clearly this is a well defined quantity which 
does not depend on the choice of the smoothing kernel ip. It remains to show that A e does indeed tend 
to A as e — > 0. To do this we let p £ {x) = X(x) — X E (x) then subtracting (28) from (23) we find that p e 
satisfies the integral equation 

Pe(x) + ^- [ D"K E (^x)p E (0^(0 = ^- [ ^K E (C,x)X(0^A0 + ^-A) [ Ao(OMO 

27T Jcr(x) 27T Jc-(x) -J JA-(x) 



2tt Jc-(x) 2tt J c -(x) 

where A is given by (30). 

Since we know from (30) that A is bounded then (24) and (25) show that 

- / \yK e {£, x)\(0 pUO + \{1 -A) [ A (Oma(0 - 0(e) 

Z7r JcT(x) 6 J A- fx) 



(31) 
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Also from Appendix C we know that C £ (x) — > C (x) as e — ► 0, and that tip — > iz r , except possibly on 
A (which is of measure zero) so that 

^ / m - [ m »no = o(e) 

*n JCo(x) ™ JC7(x) 

Hence p e (a;) satisfies the integral equation 

Pe{x) + ^-f n s K E (Z,x)p E (Z)nUO = he(x) (32) 

where h e (x) denotes the function on the right hand side of (31) and which tends to zero as e — + 0. We 
may obtain a solution to this equation by iteration and find that p e (x) also tends to zero as e — > 0. 

Thus X E (x) tends to X(x) given by (30) as e — > and inserting this into (22) and using (21) we see 
that for (p E A q for q sufficiently large 

lim(u £ ,V) = -/ (d t X(0,x a )v(x a ) + X(0,x a )w(x a ))ii a s (x a ) (33) 

Js\A 

where p,g is the volume form induced by the conical metric on S \ A and X(t, x a ) is given by (30). Thus 
U is associated to a distribution u defined by (30) and the right hand side of (33). 

It is worth remarking that even though M \ A is locally flat the Greens function for the cone is not 
sharp due to the second term in (30) (which vanishes in Minkowski space when A = 1). The solution at 
x depends not just on the initial data on Cq (x) fl S but also on points in Cq (£) fl S with £ e Cq (x) n A. 
i.e. as well as the sharp term in the Greens function there is an extra term which involves scattering off 
the axis. 

6. Conclusion 

In an earlier paper (Clarke et al 1996) we showed that it was possible to give a distributional inter- 
pretation to the curvature of a conical space-time by using Colombeau's theory of non-linear generalised 
functions. In this paper we have established that such conical singularities do not disrupt the Cauchy 
development of test fields on this fixed background. Although it is not feasible to undertake a full non- 
linear analysis and show that such singularities do not disrupt the Cauchy development of Einstein's 
equations, the higher order energy estimates obtained in Section 4 indicate that the back reaction is 
not likely to radically affect the nature of the singularity These results taken together therefore show 
that it is reasonable to interpret conical space-times not only as distributional geometries but also as 
distributional solutions of Einstein's equations. 

The concept of generalised hyperbolicity considered in this paper ((J-hyperbolicity) is different from 
that adopted by Clarke (1998) when considering curve-integrable space-times. However the curve inte- 
grate condition is primarily used to construct a geodesic congruence with tangent vector having bounded 
covariant derivative, and such a congruence may be explicitly constructed in conical space-times even 
though they fail to be curve-integrable. In a future paper we will examine the relationship between these 
two conditions. 

Finally it is worth remarking that the limiting solution that one obtains to the wave equation is more 
interesting than one might at first expect. It contains a term due to the delta- function curvature on the 
axis and therefore has tail terms even though the conical space-time (M \ A, go) is locally flat. It also 
contains a term which arises from parts of the null cone which in the limit pass through the axis and 
'mends the tear'. A naive approach in which one looked at the solution in Minkowski space in cylindrical 
polar coordinates and then rescaled the angular coordinate would not produce the correct answer. 

Appendix A. Restrictions in the Colombeau algebra 

When dealing with initial value or boundary value problems in the context of Colombeau algebras 
one is interested in solving some differential equation on R", but giving data on some lower dimensional 
subspacc S = M m . For example when solving the initial value problem for the wave equation in £(R 4 ) 

□17 = 
U\ S = V 
dtU\ S = W 
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where 5" = {(t,x,y, z) G M 4 : t = 0} = R 3 , we need to regard U G 5(R 4 ) but U| S G S(R 3 )- One way of 
doing this is to require that elements ipW G .A(R n ) are products of one dimensional kernels </? G ^4(R) so 
that 

If {/ G £(R 4 ) then we may define U ls E G(M. 3 ) by 

C/| 5 (^ 3 ) , (x, y, z)) = U(<pW , (0, y, z)) (A.l) 

Unfortunately we often want to restrict *4 g (R 4 ) to some subset which is invariant under the action of some 
symmetry group (for example rotations about the z-axis) and it is not usually possible to achieve this with 
kernels taking the form (A.l). In this appendix wc will therefore describe a way of defining restrictions 
which allows one to work with kernels with a specified symmetry. We will illustrate the approach by 
restricting £(R 4 ) to a three dimensional subspace, but the construction can be readily generalised to any 
linear subspace of K™ . 

Let / G Em (R 4 ), then we define 

j\ S : A q (R 3 ) xl 3 -»l 
by f\s(f, (x, y, z)) = /(<£, (0, x, y, z)) 

where 

ip(t,x,y,z) = \tp{x,y,z) / (<p(t, u, v) + <p(v, t,u) + <p(u, v, t)) dudv 
Jr 2 



The first important point to note is that if ip G *4(R 3 ) is invariant under some symmetry group (such as 
rotations about the z-axis) then ip G ^4(R 4 ) is invariant under the same symmetry. The second important 
feature of the construction is given by the following Lemma. 

Lemma 6. 

(1) <p G A (M. 3 ) =► <p G 

(2) cpeAqiM. 3 ) =» <peAq 

Proof. For (1) we note 



(p(t,x,y,z)dtdxdydz — | / ip(x,y, z) dx dy dz I (ip(t,u,v) + ip(v,t,u) + ip(u,v,t)) dudv dt 

= 1 

For (2) we let ip G -4g(R 3 ) with q > 1 and let a, b, c, d G N be such that l^a + b + c + ds^q. Then 

/ <p(t,x,y,z)t a x b y c z d dtdxdydz 
Jr. 4 

= | / <p(x,y, z)x b y c z d dxdy dz (ip(t,u,v) + ip(v,t,u) + p{u,v,t))t a dudv dt 

If a ^ 1, the second integral on the right vanishes, while if a = the first integral on the right vanishes. 
In either case the left hand side vanishes and (p G .4 9 (R 4 ). 

Proposition 7. 

(1) fe£ M (R 4 ) f ls e£ M (R 3 ) 

(2) / g A/"(R 4 ) =► W 3 ) 

Proof. The proof of (1) and (2) follows directly from the definition of moderate and null together with 
Lemma 6. 

Corollary 8. Let F G S(R 4 ) have representative f G £m(R 4 ), then F\s G 5(R 3 ) way 6e defined by 
F \s = [f\s]- 

We have therefore shown how it is possible to restrict elements of C7(R 4 ) to some subspace S while 
maintaining the required symmetry. The following result shows that for smooth functions this notion of 
restriction commutes with the canonical embedding i n : C°°(R") — ► Q(R n ). 
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Proposition 9. Let f e G°°(R 4 ) then 

M)\s = ts(/|s) 

Proof. 

(7|s) e = / f(0,x + ex',y + sy',z + sz')ip(x',y',z')dx'dy'dz' 
(/e)|S=g / f(et',x + ex',y + ey' , z + ez')ip(x' ,y', z') 

JKfi 

x (<p(t ; , u, v) + <p(v, t',u) + ip(u, v, t')) du dv dt' dx' dy' dz 1 
Then by looking at the Taylor expansion of (f\s) £ — {fe)\s one can verify that the difference is null. 

Appendix B. The Greens function for the wave operator □. 

In this appendix we summarise some results of Friedlander (1975) and Bruhat (1962). Let x £ M. A 
distribution may be defined on D~(x) as follows; let 

S-(x) = {tED-(x)\r(t,x)=e} 

where r(£, a;) denotes the square distance between x and £ in M. We then define a distribution 
5-(T(£,x)-e)by 



(6-(r(£,x)-e)M0)= / ¥>(0M0 

JS-(x) 

Consequently we may define a limiting distribution <5~(£, x) by 



(5-(S,x)M0) = / ¥>(0M0 

£^0+ Js~(x) 

It may be shown, for the wave operator □ (Friedlander 1975, Theorem 4.2.1), that there exists a 
function K G G°° x G°° (the biscalar) such that 

□ a;)*" (£, a:)) = z))<T(£, ar) + 2nS(t - x) (B.l) 

and moreover K can be chosen so that 

2v a rv Q A + (nr-8)A = o (b.2) 

K{x,x) = 1 

On integrating (B.l) we obtain 

/ n(K^,x)6-^,x)M0m= [ {nK(£,x))5-(£,x)tp(Q» + 2ir<p(x) 

J M JM 

which may be rearranged, by integrating by parts and using the fact that <p has compact support, to 
obtain 

V {x) = ±- [ K{£, x) L>(0 M r(0 - ^ / DK(t x) MO 

*K JC-(x) ^ JC-(x) 

This equation may be regarded as a Volterra type integral equation for ip given tp E V (Dip in this case) 

^ JC-(x) 27T Jc-(x) 

This equation naturally defines a distribution G~ : 2? — > R given by (G~,ip) — <p(x) (Bruhat, 1962). 
The distribution G~ is often known as the past Greens function. A future Greens function G+ may be 
defined by an analogous procedure, replacing the past null cone C~(x) with the future null cone C + (x) 
in all the definitions above. In fact we may define a distributions on V (g> V by 

(G ± (x,0,<p(xW(Z)) = {{E±,i > ),<p(x)) 

The fact that □ is self adjoint implies that G + (x,£) = G~(£, x). 
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Calculation of K and \3K for the wave operator □. 

We shall first consider K(x, y) in normal coordinates based at x. In normal coordinates we have 

9ab{y)y b = 9ab(x)y b 
T(x,y) = ga b (x)y a y h 



This implies that 
and so 



Substituting in (B.2) gives 
which implies 



V a r = 2y a 

□r4| fl |-V2 ao (| 5 |i/2 V - r ) 
= 2| 5 |-v^ a (| 5 |VV) 

= 8 + y a 9 a (log|.g|) 
4y a d a K + y a d a (log\g\)K = 



y a a a (M 1/4 ) = o 

This means that \g(y)\ 1 / 4 K(x, y) is constant along the geodesic connecting x and y. Therefore 

\g(y)\ 1 ^K(x,y) = \g(x)\ 1 ^K(x,x) 

= \g(x)\ 1/4 

Therefore 

K(x y) ^ ' g ^' 1/4 (B 3) 

We now consider K(x,y) in a more general coordinate system. We shall use primed coordinates to 
denote the normal coordinate system. Let 



J{y) = det 



( dy a 
\dyi» 



Then 



K(x, y) -^ /4 



ls'(y)l 1/4 

\g(x)\^J(x)- 1/2 

\g(y)\ 1/4 J(y)- 1/2 



But 



So 



~ \g(y)\ 1/4 
dr 

-Q^( x >y) = -2gab(x)y' b 

<9 2 r . , „ . ,dy' c 
:{x,y) = -2g bc {x)- 



implying that 



dy a dx b ' dy a 
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Therefore 

det(^)' 1 " 

K(x,y) = 



4|.g(.x)| \g{y)\ 

We shall now consider DK(x, y) in normal coordinates based at x. Applying □ to (B.3) 

DK=\g\- 1 / 2 d a (g ab \g\ 1 / 2 d b K) 

= \g\- 1/4 d a (g ab \g\ 1/2 d b (\g\-^)) 
= -\K(x,y)i{y) 

where 

i{y) = (g ab 9 cd 9cd, a ) b + b ab g cd g ef 9c d ,age f , b 

In normal coordinates at x it may be shown that 

g a bA x ) = 

R(x) = lg ab {x)g cd {x)g cd , ab {x) 

so 

j(x) = lR(x) 
We may now expand g a b{y) in powers of y about x; 

g a b(y) = gab(x) + \R ac bd{x)y c y d + 0(y 3 ) 

this implies 

7 (y) = lR{x) + \{^aR{x) + 2V b R ab (x))y a + 0(y 2 ) 
and on applying the Bianchi identities, 

7(2/) - lR(x) 4 0(y 2 ) 

Therefore 

UK{x,y) = -\K(x,y)R{y) + 0(y 2 ) (B.4) 
a result obtained by McLenaghan (1969). 

Appendix C: Null geodesics in a conical space-time. 

In this appendix we examine the limit of C~(x), the past null cone of x in (M, g e ) as e — > 0. To do this 
we first examine the past null cone of points in (M \ A, g ), the conical space-time with the axis removed. 
Although it is not hard to derive the equation of the null cone by working in the Minkowskian covering 
space and then making the appropriate identifications, we will derive it using a method that can also be 
applied to (M\A,g e ). 

If we work in polar coordinates the geodesics may be found using the Lagrangian 

L= 1 1 {-i 2 + r 2 + A 2 r 2 4> 2 + z 2 } 

where a dot indicates the derivative with respect to some affine parameter. The geodesic equations show 
that we may take t as an affine parameter. In view of the rotational and translational symmetries of the 
metric, without loss of generality we may consider the null cone of a point P with coordinates (0, a, 0, 0). 
We therefore consider geodesics which, when t — 0, satisfy the initial conditions: 

r = a 0=0 z = 

r = — sin 7 cos /3 (j> = sin j3 sin 7 i = cos7 
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where (3 and 7 are polar angles which parameterise the S 2 celestial sphere of null directions at P (with 
[3 = corresponding to a geodesic which hits the axis). 
Then the z equation gives 

z = t cos 7 

The null geodesies may then be found by solving the equations 

A 2 r 2 6 = J 



where 



Substituting for <p this gives 



and hence 



r 2 + A 2 r 2 d> 2 = E 2 



J = Aa sin (3 sin 7 
E 2 = sin 2 7 



dr ^sin7(r 2 - a 2 sin 2 (3) 1 / 2 
~dl ~ r 

rdr 



= -[- 

J su 



sin7(r 2 — a 2 sin 2 /J) 1 / 2 

where the minus sign is taken for the past null cone. Performing the integration and using the fact that 
r = a when t = gives 

(r 2 -a 2 sin 2 /?) 1 / 2 a cos (3 
sin 7 sin 7 

Solving for r enables one to find r{t) where 

r 2 (t) =t 2 sin 2 7-2aicos/3sin7 + a 2 (C.l) 
We now substitute for r in the <ft equation to obtain 

d(f> a sin (3 sin 7 



dt A(t 2 sin 2 7-2aicos/3shi7 + a 2 ) 

Integrating and putting in the initial condition that </> = when t = gives 

1 A sin 7 - acos3\ . 

A4>=-0±cat- 1 [ . C.2 

\ a sin /j / 

We have now obtained r(i), ^(t) and z(£) and one can now eliminate (3 and 7 and show that t, r, (f> and 
z satisfy the constraint 

a 2 +r 2 - 2ar cos(A<t>) + z 2 - t 2 = (C.3) 

The points on the past null cone of (0, a, 0, 0) are therefore the points where t < and (C.3) is satisfied. 

We now turn to considering the null cone in the regularised space-time (M, g £ ). If we use a rotationally 
symmetric smoothing kernel, as described in Appendix A, we may write the metric in the form 

ds 2 = -dt 2 + P 2 (r) dr 2 + r 2 Q 2 (r) d(f> 2 + dz 2 

where P e {r) = P(r/e), Q e {r) = Q(r/e) and P and Q are given by 

P{rf = 1(1 + A 2 ) + |(1 - A 2 )2n f (l - ^WWp 

Q(r) 2 = 1(1 + A 2 ) |(1 - A 2 )2n£(l ^)^(p)pdp 

Here <p is a smooth function with compact support such that 

2tt J (p(r)r dr = 1 

27r J Lp(r)r 3 dr = 
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so that the metric (when regarded in Cartesian coordinates as a function of ip) may be interpreted as an 
element of Colombeau's generalised function algebra £(M 4 ). 

In this space-time the geodesies may be found by considering the Lagrangian 

L £ = ±{-i 2 + P 2 r 2 +r 2 Q 2 J 2 + z 2 } 

and we see that as before we may take t as an affine parameter. We again consider the null cone of 
(0, a, 0, 0), but this time take the initial conditions: 



r = a cj> = z = 

r = — — . , sin 7 cos 6= — - - sin /3 sin 7 i = cos7 

P £ (a) aQ £ (a) 

We again have z(t) — tcos-f, and the r and t equations are now 

Q 2 (r)r 2 = J £ 
P?(r)f 2 +r 2 Q 2 (r)<j> 2 =E 2 



where 



Substituting for <j> this gives 



J £ = aQ e {a) sin (3 sin 7 
E 2 = sin 2 7 



dr_ _ ± sin 7 (Q 2 (r)r 2 ~ a 2 Q 2 (a) sin 2 /?) 1 / 2 
dt ~ P £ {r)Q £ {r)r 

and hence 



t = - 



P £ (r)rdr 



sin7(r 2 — a 2 S e (r) sin' 3 ) 1 / 2 



where S £ (r) = Q e {a)/Q £ {r). 
We now let 

R = sup{r : \tp(r)\ > 0} 

and for the moment restrict attention to geodesies for which r m ; n > eRo- Then for G A q and r > e we 
have 



e q+l 



^ = A + [^) R o) 
These estimates allow us to deduce that 

(r 2 - a 2 sin 2 /3) 1 / 2 a cos /3 , N 
sin 7 sin 7 

where g e (r) — > uniformly as e — > 0, and hence that 

r e (t) = r(t) + M*) (C.4) 

where r(f) is given by (C.l) and h E (t) — > uniformly as e — > 0. 
We now consider the equation 

d(j> _ J £ 
~dt ~ Q 2 (r)r 2 

so that 

Qe(a)cp = J 



aS 2 (r) sin j3 sin 7 



i 2 sin 7 — 2at cos /? sin 7 + a 2 + /i £ (t) 
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Using the initial conditions for <j> and fact that r £ (t) > ei? together with the estimates for P e and Q e we 
may deduce that 

4>S) = 4>(t) + k s (t) (C5) 

where <j>(t) is given by (C.2) and k e (t) — > uniformly as e — > 0. 

If we now substitute (C.4) and (C.5) into the constraint equation (C.3) for the null cone in a conical 
space-time we find we get an equation of the form 

a 2 + r 2 — 2arcos(A(f>) + z 2 — t 2 = m £ (t, r, cj>, z) 

where m e — > uniformly as e — > 0. 

We also note that as e — > 0, the excluded geodesies are those which hit the axis. So that if we denote 
by C~{x) the past null cone of x excluding those points which lie on geodesies which hit the axis, then 
the above result implies that C~ (x) — ► Cq (x) as e — > 0. 

The estimates for P e and Q e also show that the volume form /i e defined by g £ tends uniformly to a 
well defined limit fj° as e — > (at least for x ^ A). So that if € V(M) then 

hm / mm)= [ mft® (ce) 

where /if (£) and /i r (£) are the volume forms on (^"(a;) and C^x) induced by /i e and /i° respectively. 

However, as we have already observed Cq (x) is not the limit of C~(x) since it does not contain the 
part generated by geodesies which get closer to the axis than sRo, and therefore contains a tear which 
destroys the S 2 topology. We now consider the portion of the null cone generated by such geodesies 
(which restores the S 2 topology) . For a fixed value of 7 the geodesies which generate the missing piece 
are given by those which satisfy the initial condition that (3 is proportional to e. Since the geodesies 
which satisfy r m j n > sRq tend to straight lines as e — > 0, when written in quasi-Cartesian coordinates, 
the previously excluded geodesies must satisfy tan (3 < sRq /a (for sufficiently small e) . We therefore 
consider the surface S e (P) generated by past directed null geodesies n e sn (t) which emanate from the 
point P (which we may assume has coordinates (0, a, 0, 0)) making polar angles 

7G(0,7r), l3 = seR /a se[-l,l] 

and then consider the limit of this surface as e — > 0. 

Let D £ be the region for which r < eRo- Then the geodesies which generate S s leave the point 
P, enter the region D £ at Kg j7 (to) where i§ = —a/ sin7 + O(e), emerge from D e at K e s 7 (if) where 
t\ = — a/ sin 7 + 0(e), and then move outwards to infinity. By the previous analysis the two portions 
outside D e both tend to straight lines as e — > 0, while the two points «f i7 (to) and Kj 7 (tf) both tend to 
the point on the axis with Cartesian coordinates {—a/ SU17, 0, 0, —a cot 7). Thus as e — > such geodesies 
tend to a straight line with a kink in it as it passes through the axis. In Cartesian coordinates k° 7 (£) is 
therefore given by 

x = a + t sin 7 

y = 

z = t cos 7 

x = — (a + t sin 7) cos 8. 
y = — (a + t sin 7) sin 5. 
z = t cos 7 

where S s is some scattering angle which depends upon s. 

For s = 1 we have a geodesic which grazes D £ , so that its deflection will be given by that of a geodesic 
in the conical space-time, and hence S\ = —a. Similarly 5-i = a. As 7 varies between zero and ir and s 
varies between —1 and 1, then the first part of the geodesic fills in the gap in Cjj~(P) caused by geodesies 
which hit the axis, while the second part generates a portion of the surface 



> — a/ sin7 ^ t < 



> t ^ —a/ SU17 



x 2 +y 2 = (a-Vt 2 -z 2 ) 2 , t^-\z\ 
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Figure 2. 



with edges given by y = ±x tan a which exactly match with the tear in Cq (P) and restore the S 2 topology 
of the null cone. See Figure 2 for a section of the null cone through t — constant, z — constant. 

It is important to note that without imposing further conditions on the rcgularisation (such as the 
scalar curvature of g e being non-negative) we cannot assume that 5 S is a monotonic function of s or that 
it lies in the range —a ^ 5 S ^ a. However we do know that the boundary is given by Si = —a and 
5-i = a. Therefore if we consider an integral of a scalar field over the parametrised surface then (apart 
from a set of measure zero) each point outside the range occurs an even number of times (with one half 
having the opposite orientation from the other and so cancelling) and each point within the range occurs 
an odd number of times (with all but one occurrence cancelling). So that if we integrate a scalar field 
over this region we need only integrate over the region Cq (P) given by 

x 2 + y 2 = (a — \Jt 2 — z 2 ) 2 , t < — \z\, -ztana ^ y ^ stana 

which is independent of the regularisation. 

It is more natural to include the 2-surface generated by the first part of the geodesic in the definition 
of Cq(P), so that we can remove the restriction on geodesies hitting the axis and Cq(P) is now simply 
defined to be the past null cone of P in (M \ A). However (C.6) remains valid, and combining this with 
the above we obtain for all G V{M) 

1^/(^(0/4(0= />M0£r(0+ / >M0/#(0 

£ ^ ^ c o ^ c a j-j 

= / (z)v(£K(0 

Jc- 

where Cq(x) = Cq (x) U Cq (x) and /i" is the volume form induced by go on Cq (x) and Cq(x). Note 
that points on the axis are excluded from both Cq (x) and Cq (x) , but do not contribute to the integral 
on the left hand side either due to the (uniform) boundedness of fi £ even when points on the axis are 
included. 
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The final point we have to deal with is to show that A _ (x) and is well defined. A~(x) is defined to 
be the limit of C~ (x) n A as e — > 0. As we have seen k £ s 7 (ig) tends to the point (—a/ sin 7, 0, 0, —a cot 7) 
as e — > 0. So that as 7 varies between zero and 7r we see that C~(P) (~l A is given by a curve of the form 

x = 0, y = 0, t 2 - z 2 = a 2 +n s (t,z), t<0 (C.8) 

where n e — > uniformly as e — > 0. In the limit we obtain A~(x) which is given by 

x = 0, y = 0, t 2 - z 2 = a 2 , i < (C.9) 

which is independent of the regularisation. 

Also for all values of e the volume form induced on A by g e is given by the 2-dimensional Minkowskian 
value H2, (since the singular part is orthogonal to the axis and the smoothing leaves the metric parallel 
to the axis unchanged). Now by (C.9) A _ (x) is a smooth curve so the volume form induced on it 
by fj,2 is well defined and is indeed the limiting value of the volume form induced on (C.8). 
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